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Abstract. Gamma-ray emission in the TeV (10 12 eV) range has been detected by HESS from two X-ray binaries: PSR B1259- 
63 and LS 5039. In both, the early-type star provides large numbers of target photons for pair-production with TeV y-rays. 
This results in a modulation of the y-ray flux as the relative positions of the y-ray source and companion star change with 
orbital phase for the observer. The extent to which this variable absorption can provide useful diagnostics for the location 
and nature of y-ray emission is examined. The absorption spectrum and transmitted flux are calculated by integrating the yy 
cross-section along the line-of-sight, taking into account the orbit, the spectrum and the finite size of the companion star in 
LS 5039, PSR B 1259-63 and LSI +6U303, a system similar to LS 5039 but still undetected at TeV energies. In LS 5039, 
emission close to a black hole or a neutron star primary is considered. In both cases, the transmitted flux >250 GeV drops by an 
order-of-magnitude near periastron (0=0). A black hole yields a clear spectral signature in the average spectrum at ss 400 GeV. 
A neutron star yields more variability, with the spectral feature moving from 200 GeV (0=0.1) to 3 TeV (0=0.7). Only 20% 
of the flux is absorbed at 0=0.7, allowing for an almost direct view of the intrinsic spectrum. Low variability will require 
emission on large scales, more than 0.7 AU away to have <50% absorption in a jet. In LSI +61°303, significant absorption (up 
to 90% of the 100 GeV flux) is predicted only slightly before periastron, accompanied by a spectral hardening above 1 TeV. In 
PSR B 1259-63, although 40% of the flux is absorbed before periastron, the large variability seen by HESS is due to the y-ray 
emission process. The predictions made here are essential to distinguish variability in the emission of y-rays from that due to 
absorption. A modulation would provide a novel way to constrain the y-ray source. Its absence would imply that y-ray emission 
occurs on large scales. 

Key words, radiation mechanisms: non-thermal — stars: individual (LS 5039, LSI +61°303) — pulsars: individual 
(PSR B 1259-63) — gamma rays: theory — X-rays: binaries 



H 1- Introduction 



Massive X-ray binaries, composed of a neutron star or black 
hole in orbit around an early-typ e star, have lon g been recog- 
nized as possible y-ray emitters (Weekes, 1992). Recently, the 
HESS collaboration has reported very high energy (VHE) y- 
ray emission in the TeV (10 12 eV) range from two such sys- 
tems. The first, PSR B 1259-63 (=SS 2883) is a 47.7 ms ra- 
dio pulsar in a 3.4 yr orbit around a B2Ve star. Variable y- 
ray emission was detected by HESS during periastron passage 
( Ahar onian etali l2005d) . Here, the relativistic pulsar wind is 
confined by the stellar wind from the massive companion. 
Particles at the termination shock produce the non-thermal 
emission dTavani et all 1 1994 . The second, LS 5039 (=RX 
J 1826.2- 1450), is a tight 4-day binary composed of an 06V 
star and an unknown compact object. Resolved radio emission 
from the source is thought to arise from a weakly relativis- 
tic jet with v ~ 0.3c; LS 5039 has also been associated with 
the EGRET source 3 EG 1824-1514 dParedes et allkOOOl) . This 
suggests LS 5039 might be a scaled-down analogue of extra- 
galactic blazars. The VHE emission is hard (power law photon 



index a -2) an d continuously detected over a timespan of 3 
months in 2004 dAharonian et alll2005bl) . 

The large stellar luminosity (2 10 s L in LS 5039) provides 
an abundant source of seed photons for inverse Compton scat- 
tering with high energy particles. At the same time, the dense 
radiation field can absorb very hi gh energy y-ray emissio n 
through pair production yy — > e + e~ iGould & Schrederlll967l) . 
With stellar effective temperatures in the 20,000-40,000 K 
range, the typical target photon energy is a few eV which is 
in the energy range for pair production with TeV photons. The 
cross-section maximum occurs close to the threshold energy 
and has an amplitude <x yy « o"r/5 where crj is the Thomson 
cross-section. The star luminosity is very high (10 38 - 10 39 
erg s -1 ) for radii ~ 10 R Q . At a distance d ~ 0.1 AU corre- 
sponding to periastron in LS 5039 (see §3), the stellar radiation 
density in 3.5 eV photons is n* « 10 14 photons crrT 3 . A rough 
estimate of the y-ray opacity is T yy » o-y y n*d « 20. The detec- 
tion of VHE photons from LS 5039 is therefore intriguing. 

The yy opacity depends strongly on the geometry. The 
opacity can be greatly reduced compared to the simple estimate 
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above because line-of-sight y-rays travel in an anisotropic ra- 
diation field. Diffusion in th e wind is too small t o significantly 
isotropize radiation jBottcher & Dermeri 120051) . The absorp- 
tion will vary depending upon the relative location of the source 
of y-rays, the companion star and the observer. In binary sys- 
tems, the orbit and companion star type can be constrained by 
optical observations. If y-ray emission is isotropic and close to 
the compact object, absorption will be periodically modulated 
in a p redictable way JProtheroe & Stanevl Il987l iMoskalenkol 
1995). Contrast this situation with that in blazars, where the 
sources of radiation (accretion disc, broad line region) that can 
absorb y-rays close to their production site are poorly known. 
A generic overview of yy absorption in binaries is presented in 
§2. 

The effects of y-ray absorption are predicted in §3-4 for 
LS 5039, PSR B1259-63 and LSI +61°303 (= V615 Cas), a 
system very similar to LS 5039. Comparisons with existing 
or future observations (will) provide usefu l diagnostics of the 
locatio n of y-ray emission and its nature. iBottcher & DermeJ 
(2005) present calculations for LS 5039, focusing on emission 
from a black hole jet, and assuming e — and a point source 
approximation for the star. Here, expected differences in ab- 
sorption lightcurves and spectra for y-ray emission close to a 
black hole or a neutron star primary are examined for LS 5039 
and LSI +61°303 using the measured orbits and taking into ac- 
count the finite star size. Expectations if emission arises from a 
jet are considered in §5. In this work, the e + e~ pairs created by 
absorption are assumed to escape the system or radiate at much 
lower frequencies. Possible caveats when the pairs radiate back 
in y-rays and initiate an electromagnetic cascade are discussed 
in §6. 

2. VHE y-ray opacity around massive stars 

The absorption of a y-ray of energy £ on a star pho- 
ton of energy e occurs a bove an energy threshold given by 
llC,ould&Schrede 1 lfl967i) 



eE > 



2m:c* 



2„4 



(1) 



(1 -e r e*) 

where e r is a unit vector along the direction of propagation of 
the y-ray; e* is the unit vector along the direction of propaga- 
tion of the star photon (see Appendix, Fig. IA. lb . Absorption 
can only occur when the target (stellar) photons are at an angle 
to the y-ray and the minimum m^c 4 occurs for head-on colli- 
sions. 

The cross-section maximum occurs close to the threshold, 
a useful feature for quick estimates but one that errs as the an- 
gular dependence has a significant effect. For instance, y-rays 
propagating directly away from the source of target photons 
are not absorbed (the energy threshold given by Eq. ^ is in- 
finite). The decrease in the radiation field density when inte- 
grating over the photon path also reduces the total opacity. The 
opacity seen by a y-ray of given energy will therefore depend 
strongly on its emission location with respect to the star in the 
binary. 

The opacity depends also on the spectrum of target pho- 
tons. High energy y-rays are able to interact with the stellar 



Rayleigh- Jeans tail while lower energy y-rays interact with the 
exponentially decreasing Wien part. At periastron, the compact 
object in LS 5039 is only « 1 R+ away from the surface of the 
star (Fig.Q. The finite size of the star can change res ults from 
a point source approximation (e.g. Moskalenko, 1995), notably 
by giving a non-zero opacity for a radially outgoing y-ray. All 
these effects can be accurately computed (Appendix). 



2. 1 . Angular dependence 

A map of absorption around a star is shown in Fig. for 
various y-ray energies. The star has an effective temperature 
= 40,000 K and radius R+ = 10 R Q , appropriate for the 
06.5V star in LS 5039 (§3). The grey scale represents the 
amount of flux absorption exp(-r yr ) depending upon location 
of y-ray emission. The y-rays are emitted in the plane of the fig- 
ure, towards an observer situated at infinity to the right. Note 
that the solution has cylindrical symmetry about the line-of- 
sight to the star center. 

At low energies, the opacity is small until the y-ray has 
an energy above the minimum threshold for head-on interac- 
tion with stellar photons kT*. The opacity then rises quickly as 
the target radiation density initially increases exponentially. At 
progressively larger y-ray energies, the dominant effect in set- 
ting the opacity becomes the angular dependence (top panels). 

At a given y-ray energy, there is a cone pointing outwards 
close to the star for which r yy is reduced. The opening angle 
roughly corresponds to the threshold for interaction with stel- 
lar photons at the Wien cutoff. Keeping e{\ — e r e*) constant 
in Eq.H an d starting on the x-axis, y-rays interact with an ex- 
ponentially increasing density of stellar photons as the angle 
increases. However, because of the finite size of the star, the 
absorption can still be high in this cone if the y-ray is emit- 
ted close to the star. For example, >50% of the 3 TeV photons 
emitted outwards <0.5 away from the star are still absorbed 
(bottom left panel of Fig.^. 

As the energy of the y-ray further increases (bottom pan- 
els), the dominant effect in modifying the opacity is the de- 
crease in radiation density since the interaction now occurs 
with Rayleigh-Jeans photons. The decrease in opacity is then 
roughly independent of location. The situation is slightly more 
complex on-axis as the angular dependence is always important 
there. 



2.2. Orbital variability 

The existence of a cone with low y-ray opacity shows that 
VHE emission can occur close to the massive star without sig- 
nificant absorption. However, this configuration is variable as 
the binary aspect changes for the observer. Assuming isotropic 
emission close to the compact object, in a small region com- 
pared to the binary size, the resulting absorption lightcurve at a 
given energy depends on the orbital parameters (orbital period 
P or b, total mass M, eccentricity e, periastron longitude u> and 
inclination i, see Appendix). Knowledge of the orbital parame- 
ters obtained, for example, via pulsar timing or radial velocity 
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exp(-r) at 30 GeV exp(-r) at 300 GeV 




stellar radii stellar radii 

Fig. 1. Gamma-ray absorption map around a massive star. The grey-scale at each location corresponds to the integrated yy 
absorption seen by a VHE photon emitted at this location and travelling in the plane of this figure towards an observer located at 
infinity to the right. The plotted value is exp(-T rr ), on a linear scale with black corresponding to high absorption and white to low 
absorption. Absorption is shown for VHE photons with energies of 30 GeV, 300 GeV, 3 TeV and 30 TeV Full lines are contours 
corresponding to 99%, 90% (white), 50%, 10% and 1% (black) of the emitted flux being absorbed. The star is located at the 
origin, has a 10 R e radius and an effective temperature of 40,000 K. The orbit of the binary system LS 5039 is (partially) shown 
as a dotted line for illustrative purposes. Here, the binary would be seen edge-on (/=90°) by the observer at right. Gamma-rays 
emitted close to the compact object towards this observer are more or less absorbed according to orbital phase, with maximum 
absorption occurring after periastron (marked by a cross). The resulting flux modulation would be close to that shown in Fig. [5] 
for the case of /=60°. 



studies, leads to constraints on the physics of VHE emission by 
comparing expected and observed lightcurves. 

The case for which absorption is invariant of orbital phase 
is for a face-on, circular (i = e = 0). An inclined, circular orbit 
results in absorbed flux lightcurves which are symmetric with 
respect to the orbital phase of superior or inferior conjuction. 
The changing orbital separation in an eccentric orbit varies the 
absorption even with i — 0°. Eccentric orbits typically result in 
y-ray lightcurves with a flux dip at periastron. 

2.3. Absorbed spectrum 

The maps in Fig. show that the extinction depends on 
location but also on energy. Pure y-ray absorption pro- 
duces an absorption trough in the spectrum that, if mea- 
sured, can provide important di agnostics for emission pro- 
cesses ( Bottc her & Dermeri 120051) . The location of maximum 
absorption changes as the compact object revolves around the 
massive star, providing further diagnostic power. The trough is 



quite wide, becoming obvious only for data covering a large 
span in energies. For temperatures typical of massive stars, this 
span is ~ 30 GeV-30 TeV (§3). The emitted y-ray spectrum is 
softened by absorption in data limited to the low energy end 
(<300 GeV) while it is hardened at high energies (>300 GeV). 
In contrast, the spectral energy distribution of the extragalac- 
tic background light results in a softeni ng at TeV energies o f 
the emitted y-ray spectrum from blazars (Steck er at all fl992>. 
Gamma-ray absorption by Galactic background radiation is 
negligible for X-ray binaries. 

3. Gamma-ray absorption in LS 5039 

The origin of the VHE emission detected by HESS from 
LS 5039 has yet to be established. An intriguing possibility 
is that the emission occurs in a relativistic jet as in blazars. 
Resolved radio emission in LS 5039, attributed to a mildly 
relativistic jet, constitutes a strong argument in favour of 
this possibility. Alternatively, the VHE emission could be a 
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result of the interaction of a pulsar relativistic wind with 
the stell ar wind from t he massive star, as in PS R B1259- 
63 jMaraschi & Treves! Il98ll iTavani et all Il994 . The sim- 
ilarity in spectral energy distributions between LS 5039 and 
PSR B 1259-63 argues in favour of this interpretation. 

Can the expected y-ray absorption be used to investigate 
the origin of the y-ray emission ? The star temperature and 
radius derived fro m optical spectroscopy are * 39, 000 K 
and R* * 9.3 R G JCasares et alJl2005al thereafter C05). The 
star mass is s* 23 M . Radial velocity studies provide a good 
determination of f or b ~ 3.906 days, eccentricity e « 0.35 
and periastron angle a> » 226° . The orbital so lution of C05, 
which includes data from M cSwain et alJ J2004I) whilst adding 
better phase coverage, supersedes previous determinations. By 
assuming that the secondary is pseudo-synchronized (in co- 
rotation at periastron), C05 find i m 25° and a compact ob- 
ject mass 3.7 + 1.3 M . Whether the star has had time to 
pseudo-synchronize in this young, massive X-ray binary is not 
clear, and a neutron star cannot be firmly ruled out yet. Taking 
into account the mass function and the more conservative con- 
straints on the system inclination i derived by C05, two orbital 
solutions are adopted here to illustrate the possible extremes: 
one with i = 60°, corresponding to a canonical 1.4 M neutron 
star; and one with i - 20°, implying a 4.5 M black hole. 

Hereafter, the y-ray emission is assumed constant, isotropic 
and occuring in a small region (compared to the binary separa- 
tion) close to the compact object. Pairs created by the absorp- 
tion of y-rays are further assumed to escape freely or to radiate 
their energy at frequencies below those of interest. These as- 
sumptions are examined in §5-6. The aims are to investigate 
(1) how much absorption occurs and at what phases; (2) ex- 
pected observational signatures; (3) whether a black hole or a 
neutron star primary can be distinguished. 

3.1. Orbital variation of absorption 

The orbital variation of y-ray absorption at various energies is 
shown in Fig.[5]for i - 20° and i = 60°. In both cases, signifi- 
cant absorption in the 0. 1-1 TeV band is predicted at periastron. 
The opacity reaches r yy « 40 at phase <p j» 0.1 for 200 GeV 
photons (z = 60°). However, the opacity is smaller than unity 
over a large fraction of the orbit. VHE y-rays escape because 
of the anisotropy of the target radiation field. The minimum ab- 
sorption (peak flux) occurs at (f> « 0.6 when the compact object 
is far from the star and close to the observer. 

A non-detection of the periodic variation by HESS would 
suggest that y-ray emission occurs far from the companion star. 
Indeed, the available HESS lightcurv e from 2004 is statistically 
constant at the 15% confidence level ( Aharonia n et allE0 05b). 
A pe riodogram s hows n o peak with more than 2<x significance 
( Aharoni an et all l2005bh . However, only 14 nightly flux aver- 
ages over 120 days are available. Future observations will allow 
firmer conclusions to be drawn. 

The amplitude and asymetry of the variation decrease with 
inclination. However, these differences may be difficult to dis- 
tinguish when comparing to observations of lightcurves inte- 
grated above the instrumental enegy threshold. The absorbed 
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Fig. 2. Modulation of the y-ray transmission exp(-r ry ) with or- 
bital phase in LS 5039 for two values of the system inclination 
i. A canonical neutron star corresponds to i = 60° while i = 20° 
implies a 4.5 M black hole. Isotropic y-ray emission is as- 
sumed to occur at the location of the compact object. The solid 
line is for 100 GeV photons (H.E.S.S. threshold), dashed is for 
1 TeV and dotted for 10 TeV photons. Periastron passage is at 
<* = 0. 



fluxes for the two inclinations a re plotted in Fig. [3] along with 
the HESS 2004 lightcurve from I Aharonian et alJ iEo05bl) . The 
intrinsic spectrum was assumed to be a power law dN y oc 
E~ 2A dE, which was integrated after absorption above the same 
threshold as the HESS data (250 GeV) to yield the expected 
flux variations. The average expected photon index is consis- 
tent in both cases with the HESS measurement (-2.12+0.15, 
lAharonian et al.l2005bl) . For an inclination of 20° (black hole), 
the average fraction of the intrinsic flux that is absorbed above 
250 GeV is 89%. For i = 60° (neutron star), the fraction 
is 64%. In Fig. [5] the observed and expected lightcurves are 
renormalised to their maximum. The length of the absorption 
'eclipse' is longer and the peak flux occurs later in phase for a 
high inclination. 
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Fig. 3. HESS o bserved lightcurve of LS 5039 in 2004 (points 
with error bars, lAharonian et alJEo05bl) . folded on the orbital 
period using the ephemeris of C05. The expected lightcurve 
caused by absorption is plotted as a solid line for i = 60° and 
as a dashed line for i = 20°. Here, the lightcurves are integrated 
above 250 GeV and normalized to their maximum (Fig. |2 
showed the absorption for given energy). The orbital phase of 
the peak and the duration of strong absorption (> 90%) are dif- 
ferent between emission close to a black hole (i = 20°) and 
close to a neutron star (i = 60°). 

3.2. Absorption spectrum 

The average absorption is greater for low inclinations as the 
angle between the star, emission location and observer (t^, 
Fig. IA. li stays close to 90°. Therefore, spectral signatures of 
yy absorption should be easier to see in a phase-average spec- 
trum if emission occurs close to a black hole. 

The absorption spectra expected at various orbital phases, 
and for the two inclinations considered, are shown in Fig. |4] 
There is a clear absorption at all orbital phase with limited 
spectral changes in the i = 20° case (black hole). The shape 
is roughly independent of phase with maximum absorption oc- 
curing around 400 GeV. In contrast, the energy of maximum 
absorption changes dramatically in the i = 60° case, shifting 
from 0.2-0.3 TeV to 2-3 TeV when going from orbital phase 
<p = 0.1 to <p — 0.7. Hence, more dramatic changes are expected 
in phase-binned spectra if the inclination is high. 

The 2004 HESS spectrum covers the 0.3-3 TeV range. 
Weighting by the integrated flux, the phase-averaged absorp- 
tion for i = 20° would produce a hardening of almost 0.5 of the 
intrinsic photon index between 0.3-3 TeV, assuming power-law 
emission at the source. Of course, a single power-law fit would 
not be an adequate description for a high S/N spectrum of the 
absorbed y-rays. When i = 60°, the amplitude of the variations 
in absorption is greater, but the phase-averaged hardening is 
slight, » 0. 1 in photon index. Spectral changes are complex in 
this case, depending upon energy threshold. Above a few TeV, 
the overall effect is a spectral hardening at low fluxes compared 
to high fluxes. With a threshold at a few 0. 1 TeV, the spectrum 
may appear to soften at high fluxes as the absorption trough 
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Fig. 4. Spectrum of y-ray absorption at selected orbital phases 
<p for two inclinations i of LS 5039 (see Fig. The absorp- 
tion spectrum varies little for i = 20°. The line-like absorption 
feature may therefore be easier to detect for low inclinations, 
despite the stronger absorption near periastron for i = 60°. In 
the latter case, although spectral variations will be easier to de- 
tect, the large differences in absorption spectrum can swamp 
out the line feature in a phase-averaged spectrum. An almost 
absorption-free view of the intrinsic spectrum is expected at 
flux maximum for i = 60°. 

moves to more than a TeV and more low energy photons are 
detected. For a neutron star, the orbital phase at which the flux 
peaks provides an almost direct view of the intrinsic spectrum, 
with at most «20% absorption at 2-3 TeV. The intrinsic spec- 
trum is always significantly absorbed (r rr > 1) in the black 
hole case. 

4. Gamma-ray absorption in LSI +61 303 and 
PSR B1 259-63 

Absorption in LSI +6T303 and PSR B1259-63 is investi- 
gated here under identical assumptions as LS 5039 for the 
gamma-ray emission (constant, isotropic, close to compact ob- 
ject). PSR B1259-63 has been detected by HESS near peri- 
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astron. Inaccessible to HESS (but accessible to MAGIC and 
VERITAS), LSI +61°303 has yet to be detected by ground- 
based Cherenkov telescopes. 



4.1. LSI +61° 303 

LSI +61° 303 is a massive X-ray binary with a very simi- 
lar spectral energy distribution to LS 5039. It also has re- 
solved radio emission. LSI +61°303 has long been associated 
to hard y-ray emission with possi ble counterparts seen by C os- 
B and EGRET. This prompted iMaraschi & Treves! Jl98ll) to 
propose that high energy emission in LSI +61°303 is powered 
by an interacting pulsar wind. Jet m ode ls have been propose d 
bv lBosch-Ramon & Paredesl d2004 and lRomero et alJ d2005t) . 
The BOV star has T+ * 22, 500 K and fl* * 10 Rq. Recently, 
ICasares et alJ J2005bl) have updated the orbital parameters from 
optical spectroscopy, finding P or b ~ 26.496 days, e * 0.72, 
to ss 21°. Radial velocity studies are insufficient to rule out a 
neutron star. As with LS 5039, two cases are studied, represent- 
ing extremes compatible with observations; one with i = 60° 
for a 1.4 M© neutron star, and one with i = 20° for a 4 M Q 
black hole. 

The absorption lightcurve is shown for both cases in Fig. [5] 
The orbit is wider than in LS 5039 so that the variation in y- 
rays absorption occurs essentially only close to periastron. The 
distance to the star is then w 0.1 AU, as with LS 5039. The 
stellar luminosity is lower in LSI +6T303 than in LS 5039, 
so that not all the y-rays are absorbed at periastron despite the 
similar distance. At apastron the orbital separation is « 0.7 AU 
but there is still some absorption, around 20% of the 1 Te V flux, 
because the compact object is behind the star when projecting 
the orbit onto the sky. 

The amplitude of the variations is greater at high inclina- 
tions, as in LS 5039. The spectral changes expected close to 
periastron for i = 60° are dramatic, with peak emission that fol- 
lows the minimum within a few days. Fig. [5] shows the strong 
hardening of the spectrum above 1 TeV with lower fluxes. 
The absorption feature outside of periastron is a 20% effect. 
LSI +61°303 suffers less absorption for most of its orbit so 
should be more easily detected at VHE energies than LS 5039, 
if the same processes are indeed at work in the two. 



4.2. PSRB1259-63 

Pulsar timing in PSR B 1259-63 allows an accurate determi- 
nation of orbital parameters: P or b ~ 1236.7 days, e « 0.87, 
(star) to * 138.7°, i « 35°. The star has « 27, 000 K and 
R+ » 10 Rq for a total luminosity of 5.8 10 4 L (Lfohnston et all 
Il994l) . Fig. Q shows the expected lightcurve around perias- 
tron passage, taking only absorption into account. The y-rays 
are emitted at the pulsar location with a power law spectrum 
dN y oc E~ 2 *dE. The pulsar is w 10 13 cm at its closest to the 
Be star, further away than LS 5039 or LSI +61°303. Although 
the pulsar is >20 R+ away from the star, there is still some ab- 
sorption because of the orbital orientation; photons emitted at 
<p x must travel near the star on their way to the observer. 
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Fig. 5. Modulation of the y-ray transmission with orbital phase 
in LSI +61°303 for i = 20° (4 M black hole) and i = 60° 
(1.4 M neutron star). The solid line is for 100 GeV photons, 
dashed is for 1 TeV and dotted for 10 TeV photons. Significant 
absorption occurs only slightly before periastron (<p = 0) in 
both cases. As with LS 5039, a higher inclination leads to a 
larger amplitude of absorption. 



Maximum absorption occurs slightly before periastron, 
with up to 40% of the y-rays above the HESS threshold 
(380 GeV) being absorbed. The measured HESS lightcurve has 
a greater amplitude with a minimum occuring later than ex- 
pected. Absorption is clearly not sufficient to explain the obser- 
vations and some of the variability must be intrinsic. A likely 
possibility is the angle-dependence in the inverse Compton in- 
teractions producing the TeV emission ( Kir k et allll999ft . VHE 
emission would not be isotropic. Models should be compared 
to absorption corrected fluxes. Finally, hardening by absorption 
of the 0.3-3 TeV spectrum is at most 0.2 in photon index, too 
small to be measured by HESS. 

5. LS 5039 and LSI +61°303: opacity in a jet ? 

Are there good reasons to expect that y-rays are emitted 
at the location of the compact object? For a binary plerion 
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Fig. 6. LSI +61°303 absorption spectrum at orbital phases 
4> = 0.05,0.5 and 0.95 assuming i = 60° (see bottom panel 
of Fig.|5jl- There is a dramatic spectral change within a couple 
of days of periastron. 
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Fig. 7. Modulation caused by y-ray absorption around perias- 
tron in PSR B 1259-63. The data points are the HESS integrated 
lightcurv e above 380 GeV measu red during the 2004 periastron 
passage ( Aharonia n et alll200 5c). The data are renormalised to 
the maximum 10~ 7 irT 2 s _1 . The solid line is the integrated flux 
above the same threshold after absorption, assuming constant 
power law emission wih a photon index of -2.8. Absorption 
plays a minor role in the observed variability. 

like PSR B1259-63, most of the VHE em ission should occur 
close to the pulsa r wind termination shock ( Tav ani et all 1 1994: 
iKirk et all fl99^ . The overall geometry is set by the stagna- 
tion point R s , the point on the orbital axis wher e ram pressure 
from the stellar and pulsar wind balance ( Mar aschi & Trevesl 
Il98lh . For LS 5039, (R s /d) ~ O.l^f where d is the orbital 
separation (0.1-0.2 AU in LS 5039) and E 36 = 10 36 erg s 1 
is the pulsar spindown power (taken equal to that measured in 
PSR B 1259-63). The assumed stellar wind has a mass loss rate 
of 10~ 7 M yr~' and a coasting speed v w = 2000 km s _! , appro- 
priate to LS 5039 (C05). The stagnation point in PSR B 1259-63 



is als o at about a tenth of the orbital separation dTavani et all 
1994). VHE emission in this case can be expected to occur on 
a small scale compared to the orbital separation, close to the 
pulsar. 

There is no such certainty with a black hole (or neutron 
star) jet. Jet models for the y-ray emis sion have been prop osed 
for both LS 5039 and LSI +61°303 dBosch-Ramon & Paredest 
l2004lRomero et all l2005UParedes et alll2005l) . Although the 
gravitational potential energy reservoir that can be tapped is 
largest close to the black hole, most of the energy may not 
be radiated there if it is carried by protons or Poynting flux 
(Sik ora et all 120051) . Shocks or reconnection can convert ki- 
netic or electromagnetic energy into radiation progressively 
over large distances, up to a termination shock with the ISM. 
There is evidence for X-ray synchrotron emission on parsec 
scales from high ener gy electrons in the relativistic jet of 
XTE Jl 550-564 dCorhel et alll2002l) . 

5.1. Variation with height 

The change in flux absorption with height above the orbital 
plane in LS 5039 is shown in Fig. [8] The inclination is 20° 
and the black hole is at periastron. Absorption decreases to 
less than 50% at all energies when emission occurs more than 
» 0.7 AU awa y from the plane, as d i scusse d in lAharonian et all 
d2005albl) and iBottcher & Dermerl J2005h . The same calcula- 
tion for LSI +61°303 yields a distance of « 0.1 AU. If VHE 
emission can be independently linked to a relativistic jet, a lack 
of temporal or spectral absorption signatures will imply emis- 
sion at large distances. 

Unlike TeV blazars, the target photon field here is external 
to the jet and yy absorption will not be reduced by relativistic 
aberration. At a given energy, the dependence of the opacity 
with height is se t by the c ompanio n light, prov iding constraints 
on models as in lBlandford & Levinsonld 19951) where the jet is 
structured by the y-photosphere. 

5.2. Beaming 

The height at which r yy < 1 in LS 5039 corresponds to roughly 
10 7 gravitational radii R g for a stellar mass black hole. In 
blazars, VHE emission is often proposed to occur in the broad 
line region at parsec-scales, i.e. about 10 4 R g for a 10 9 M black 
hole. In the internal shock model, the distance of first interac- 
tion is ~ Y 2 cAt where At is the observed timescale a nd T i s 
the typical bulk Lorentz factor of ejecta dSpada etall 1200 ll) . 
Assuming LS 5039 is a genuine microblazar with scaled-down 
physics, then At would be reduced with black hole mass. For 
the given emission distances, the conclusion would be that F in 
LS 5039 is at least the same, if not higher than in blazars. 

Could ejections in LS 5039 have T » 10, as in- 
ferred in blazars from opacity arguments and y-ray modelling 
dDondi & GhiselliniL Il995l iGhisellini et ail 1 19981) ? Emission 
from a jet will be c/eboosted if viewed at an angle larger than 
« 25° for F=10. The system inclination required for a black 
hole primary is close to this value. In this case, emission from 
a jet perpendicular to the orbital plane would not be boosted. 
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Fig. 8. Variation in transmitted flux with height z above the 
plane in LS 5039. The assumed jet is perpendicular to the or- 
bital plane, the inclination is 20° and the black hole is at peri- 
astron. The absorption is shown for 100 GeV and 10 TeV pho- 
tons. The curve for 1 TeV photons would be indistinguishable 
from the 100 GeV curve in this plot. VHE photons emitted a 
few AU away from the black hole suffer little absorption. For 
comparison purposes, the orbital separation in LS 5039 varies 
between 0.1 and 0.2 AU. 

This might fit with the weak X-ray to y-ray luminosity ob- 
served (~ 10 34 erg s _1 ). Then, LS 5039 would not be a genuine 
microblazar with the jet along the line-of-sight. 

Current VLBI radio observations of LS 5 039 only indi- 
cate a mildly relativistic outflow with v < 0.4c ( Pa redes et all 
2000). Yet, the Lorentz factors of TeV blazars deduced from yy 
absorption are known to be much greater th an those measured 
at lar ger distances by VLBI observations dPiner & Edwards! 
120041) . In both cases, there would be a puzzling discrepancy 
between the high F inferred from y-ray studies and the low T 
measured on larger scales in radio. 

6. Effect of pair production 

The e + e~ pairs created by absorption can boost stellar pho- 
tons to very high energies via inverse Compton interactions. 
These y-rays can in turn be absorbed and a cascade follows. 
The initial pair created by absorption of a 1 TeV photon has a 
Lorentz factor y « 10 6 . In LS 5039, their inverse Compton 
interaction with stellar photons occurs in the Klein-Nishina 
regime on a timescale ?kn ~ 8y6^ [ [In y& + 1 .2] seconds 
(Blu menthal & Gouldl fl970). with d .i=0.1 AU. Hence, the 
typical distance travelled by the pair before a photon is boosted 
to high energies is small enough compared to the binary size 
to start a cascade. However, the cascade will be quenched if 
the pair loses energy preferentially to another emissi on process 
that does not re-emit VHE y- rays (e.g. synchrotron: iBednarekl 
ll997HAharonian et all2005al) . 

The cascade reduces the effective r yy , since absorbed y- 
rays are partly re-radiated (at lower energies). If T yy is very 
high, energy initially radiated in the absorbed band is de- 



graded to frequencies below the pair production threshold (< 
30 GeV in LS 5039). The pile up at these frequencies might 
account for the luminosity step between the EGRET (GeV en- 
ergies) an d HESS (TeV energies) measurements (see Fig. 1 in 
lAharonian et al.lEo 05a). Absorption of a constant spectrum in 
vF v from GeV to TeV energies cannot account for this differ- 
ence. Such a spectrum would, on average, only be absorbed 
by a factor 3-10 above 250 GeV (§3.1) when the luminos- 
ity jump is »30. Taking into account the photons degraded to 
lower energies by the cascade would raise the absorption step. 
Orbital variability would then even be expected at energies be- 
low the interaction threshold (< 30 GeV, hence in the EGRET 
or GLAST range ), in antiphase with the absorption lightcurves. 
Bednarek ( 1997) called this "y-ray focusing by soft radiation". 

Angular effects must be considered for yy opacities of a 
few. The estimate of the Klein-Nishina interaction timescale 
assumes an isotropic photon bath. Just as the yy absorption 
is changed by the geometry, the pair motion with respect to 
the star will change the inverse Compton opacity. Maximum 
boost will be achieved for head-on collisions with stellar pho- 
tons. Because of momentum conservation, the cascade should 
tend to develop along the direction where the initial photon 
was emitted. However, the cascade is effectively isotropized 
as the gyroradius is small compared to the interaction length 
dBednarekUl997l) . Radiation can then be redistributed at other 
orbital phases. 

A quantitative description o f these e f fects is fea- 
sible by Monte-Carlo meth ods dBednarekl 1 19971 [2000; 
Sie roowska & Bednarek! l2005h but is beyond the scope of this 
work. A model of the intrinsic spectrum should be considered 
before doing this. Cascades will probably be important only 
when the yy opacity is high and the source behind the star in 
the plane of the sky. This ensures maximum energy boost and 
maximum density of stellar photons along the initial develop- 
ment of the cascade. This is the case at periastron in LS 5039 
but not in LSI +61°303 (source close to front of the star at max- 
imum absorption) and even less in PSR B1259-63 (low r yy ). 
The amplitude of flux variations will be smaller, though per- 
haps not noticeably so on a linear scale and with energy thresh- 
olds > 100 GeV, which are usually at or above the energy of 
maximum absorption (e.g. compare the Ce n X-3 lig htcurves 
above 300 GeV with and without cascade in Bednarek 2000). 
The absorption spectrum around periastron will be softened. 
The changes will be small when r yy < 1. Qualitative differ- 
ences in absorption at TeV energies for a neutron star or black 
hole case should therefore remain. 

7. Conclusion 

The intense stellar radiation field in massive X-ray binaries pro- 
vides a novel way of investigating very high energy processes 
around compact objects. Because of the dependence of the yy 
cross-section on angle, the anisotropy of the field has a strong 
effect on the opacity to very high energy y-rays > 10 GeV. The 
orbital orientation to the observer and location of emission are 
crucial in determining the absorption. In PSR B 1259-63 this 
leads to significant absorption before periastron, although the 
pulsar is > 20R* from its stellar companion. 
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The changing radiation field modulates the gamma-rays in 
a predictable fashion. The result is a controlled, repeatable 
experiment that can help pinpoint the location of emission. 
Comparing the predictions to observations is essential to sin- 
gle out the variability at the source from that due to absorption 
in the binary. In LS 5039, the binary separation is only 2-4 /?* 
yet absorption can be much smaller than implied by simple es- 
timates. Indeed, absorption of y-rays emitted close to a neutron 
star (implying a high inclination) is negligible at inferior con- 
junction. The spectrum above 1 TeV is considerably hardened 
near periastron in LS 5039 and LSI +61°303. LSI +61°303 is 
largely absorption-free which should ease its detection at high 
energies by MAGIC or VERITAS, if similar processes are at 
work in both sources. 

The detection of an orbital modulation of the TeV flux in 
LS 5039 would imply emission within the system, either from 
a pulsar wind termination shock or from the inner region of a 
jet. If the primary is a black hole, the inclination will be » 20° 
and the absorption spectrum has a break at * 400 GeV with 
limited orbital spectral changes. For a neutron star, the incli- 
nation is w 60° and the break is expected to change from 0.2 
to 2 TeV from low to high fluxes. The absence of orbital vari- 
ations or spectral signatures in VHE observations of LS 5039 
would suggest emission occurring in regions of low yy opac- 
ity. Jet emission at large distances compared to the system scale 
would be the prime candidate. 
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Appendix A: Calculation of y-ray absorption in a 
massive X-ray binary system 

A. 1 . Opacity in vicinity of the star 

The differential absorption opacity seen by a y-ray of en- 
ergy E located at a position P and travelling in the direc- 
tion e r due to photons o f energy e emitted at S along e* is 
JGould & Schredeilll967l) 

dTyy = (1 - e y e*)n E cr yy dedQd/ (A.l) 

e y and e* are unit vectors in the directions of propagation; dQ 
is the solid angle of the surface emitting the e* photons; n £ is 
the radiation density in photons cirT 3 erg -1 sr _1 . 

The yy interaction cross-sectio n <x yy depends only on s, de- 
fined as JGould & Schrederlll967b 

S= 2^ (1 - e ^ (A ' 2) 

and the photons interact only if s > 1. The variable s can be 
rewritten as s=e/e m ; n where e m ; n is the threshold for interaction 
at y-ray energy E (Eq.Q. 

The massive star dominates the radiative output in LS 5039, 
LSI +61°303 and PSR B1259-63 so that any other source of ra- 
diation for pair-production with TeV y-rays is neglected. The 
star has a radius and is assumed to have a blackbody radia- 
tion density with a temperature 

^ = f? exp (e /^)-l Ph ^ 3erg " S1 '" (A3) 
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Fig.A.l. Geometry for absorption of a y-ray at a location P 
due to e + e" pair creation with stellar photons emitted at S . The 
y-ray is emitted at / and / is the length of the y-ray path to P. 
The (x,z) plane is defined by the star center and y-ray path. 

Orienting the coordinate system so that the y-ray path is in 
the (x,z) plane ("Fig. IA.lt . then the location and direction of the 
y-ray are given by the angle if/ and distance d to the star center 
and 



1 - e y ■ e* = 1 + cos if/ cos 9 + sin if/ cos <p sin 9 



(A.4) 



The solid angle of the star surface element emitting photons 
along e* is 



dQ = sin0d0d6» 



(A.5) 



The length / of the y-ray path since its emission at a distance 
do and at an angle if/g (point I in Fig. lA.ll with if/o-Q for a y-ray 
travelling directly towards the center) is related to if/ so that: 



if/ = tan 



1 do sin t//o 
ydo cos if/o - I 



for / < do cos if/o 



(A.6) 



and if/ = 7r + tan l (...)forZ > docosif/o- The distance d is related 
to I by 



d 2 = d\ + I 2 - Idol cos iJ/q 



(A.7) 



The dependence of the differencial opacity dr yy on E, e, 
4>, 9 and / is now explicit. The total opacity to infinity for a y- 
ray of given energy E and initial location requires a quadruple 
integral: 



' 77 



I dl I dcosfl ( dip ( 

J0 Jc min Jo Je„ 



dr 



yy 



dedQdZ 



;de 



(A.8) 



where c m ; n = (1 - R 2 /d 2 ) 1 ^ 2 . It is convenient to rewrite the en- 
ergy integral as a definite integral on/3 = (1 - 1 /s) 1/2 within the 
interval [0, 1]. Similarly, the integral on I can be transformed 




1.0 1.5 2.0 2.5 3.0 

distance (R.) 

Fig. A.2. For a y-ray emitted directly away from the star (t/r = 
n), the point source approximation (Eq. IA.9> predicts no ab- 
sorption. The exact transmitted flux exp(-r yy ), calculated us- 
ing Ea. IA.8l is shown here for a 1 TeV y-ray as a function of 
the distance of its emission point (/?*=10 R , T* =40,000 K). 
Because of the star's finite size, there is absorption even for 
tfr — n when the y-ray is at distances < 2R*. 



into a definite integral on if/ between [if/o,n] with the condition 
do sin if/Q > R+ (or the y-ray hits the star). An iterated Gaussian 
quadrature is used to calculate the integral. 

When R+ <sc d, an integration of Ea. lA.8l over 9 and <p re- 
turns the point source approximation 



cr yr « e (l + cos if/)dlde 



(A.9) 



The threshold is e m ; n = 2m 2 c 4 /E(l + cosif/). Numerical com- 
parisons show the approximation works reasonably well except 
for if/ x n (Fig. IA.2t or when d - R+ < R+ at closest approach. 
The full calculation is then preferable. The results presented 
here were derived using the exact calculation. 

A.2. Orbital dependence 

Assuming emission close to the compact object, the initial dis- 
tance do of the y-ray to the star is simply 



do 



a(l-e 2 ) 



1 + e cos(# - (jS) 



(A. 10) 



with the semi-major axis a = (GMP 2 rh /4n 2 ) 1 ^ 3 (M is the to- 
tal mass), e the eccentricity, 9 the true anomaly and a> the 
periastron angle. The orbital phase (mean anomaly) is <p — 
(rj — e sin T])f27T where 77 is the eccentric anomaly with 



77 1 - e 

tan - = 

2 \l+e 



1/2 



tan ■ 



(A.ll) 



The angle if/o of emission to the observer (located at 00) mea- 
sured from the star, at each point of the orbit (Fig. IA.lt is 



cos if/o = sin 0sin i 



(A. 12) 
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where i is the inclination of the orbit. If emission occurs at a 
height z above the location of the compact object and perpen- 
dicular to the orbital plane (as in Fig. [3), then the initial distance 
to the star simply becomes d = (d^ + z 2 ) 1 ^ 2 and the initial an- 
gle is changed to cosiA = (do /d) cos iff (with do and if/ given 
above). 



